Phase diffusion of a two-component Bose-Einstein condensates: exact and short-time 

solutions for arbitrary coherent spin state 
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We investigate phase diffusion of a two-component Bose-Einstein condensates prepared initially in 
arbitrary coherent spin state \6o, <f>o). Analytical expression of the phase-diffusion time is presented 
for 8o 7^ 7r/2 case. In comparison with the symmetrical case (i.e., 9q = 7r/2), we find that the 
diffusion process becomes slowly due to the reduced atom number variance. 
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Phase diffusion of Bose-Einstein condensates (BECs) 
043 1 aroused from atom-atom interactions destroys 
phase coherence, and thus sets a limit to the applica- 
tions of the condensates in high-precision measurement 
and quantum information processing. The phase diffu- 
sion is observable by measuring fringe visibility in atomic 
interference experiments [8Hl2j. For a two-component 
BECs prepared initially in a coherent spin state (CSS) 
\e ,<p ) with O = tt/2 [see Ref. [H, or Eq. ©], it 
has been shown that the single-particle coherence (or 
fringe visibility) decays exponentially with the time scale 
Xtd = (N/2)~ 1 ' 2 [2j-|5(, dependent upon atom number TV 
and the self-interaction strength x ~ ( a n+ a 22 — 2ai2)/2. 
Here, an and a.y (i ^ j) denote intra- and inter-species 
s-wave scattering lengths, respectively. In this paper, we 
generalize previous works by considering arbitrary initial 
CSS, i.e., 8o 7^ 7r/2. Exact and short-time solutions of 
the single-particle coherence are obtained. In comparison 
with 9q = 7r/2 case, we find that the diffusion becomes 
slowly due to relatively narrow atom number variance. 

Considering a two-component BECs with internal 
states |1) and |2) confined in a deep potential, we adopt 
single- mode approximation (SMA) [1J], i.e., keeping the 
condensed- mode wave function $o( r ) for the two com- 
ponents, so the total system can be well described by 
second-quantized Hamiltonian (K = 1): 



H 



U)2ll 2 ~\ 
2/1 ^2 



(i) 



where hi, b\, and hi (= h\bi) are the annihilation, cre- 
ation, and number operators for the internal states i = 
1,2, respectively. In addition, Wj are single-particle ener- 
gies, and Uij = (Airaij/M) J d 3 r|$ ( r )| 4 are atom-atom 
interaction strengthes. Particle number operator A" = 
hi +ri2 is a conserved quantity and is set to the c number 
N. Hamiltonian ([1]) can be rewritten as H — 5J Z + xJ%i 
where the detuning S — u; 2 — uji + (u22 — un)(N— 1)/2, and 
the self-interaction strength x = (1*11+^22 — 2u\2)/2, The 
spin operators J+ = (J-)^ = b\b\ and J z = (h2 — h%)/2, 
obeying SU(2) algebra. The nonlinearitvY J? has been 
proposed to prepare spin squeezed state 15l-ll7| and also 
quantum entangled state 3 1!|. Note that the model 



considered here with nonzero S can be also used to study 
the BECs in an asymmetric double- well [2(J. Squeezing 
via coupling of the BECs in a double- well potential with 
a cavity light field has been investigated in Ref. 
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Besides the spin squeezing, mean-field interaction X-^t 
also leads to phase diffusion in the two-component BECs 
[2-0, [ijj], which can be illuminated schematically by 
Husimi Q function Q(6,<f>;t) = \ (8, </>|*(i)) | 2 , where 

\6,4>) = exp{ j6»(J x sin </)- J y cos</>)}|j, j) (2) 



is arbitrary CSS [1J] and |*(t)) 



-iHt 



*(0)) is a state 



vector at any time t. For an initial state 1^(0 
analytic expression of the Q function reads ll 



Q(M;0) = \(e,0\e o ^o)f 



1 + cos e 



2.; 



(3) 



where cos 9 = cos 8 cos 8q + sin 8 sin 8q cos(</> — 4>q). The 
Q function can be plotted in three-dimensional phase 
space (i.e., Bloch sphere) [ijj Il6j . or alternatively, in 
a two-dimensional phase space (<p,s z ) via a mapping 
Sz = jcosd [ll EE HI- As shown in Fig. []]>) and 
Fig. Qlb), the density of Q{8, <fr; 0) is distributed isotrop- 
ically, indicating the minimal uncertainty relationship of 
the initial CSS. Under the government of nonlinear inter- 
action yJ|, the system will evolve into a spin squeezed 
state 15] with anisotropic distribution of the Q function 
[see Fig. QTcVl . As time increases, the spin system is over- 
squeezed [HI, and shows a spread of the Q function 
along the <p axis, which simulates an increased relative 
phase fluctuation (i.e., phase diffusion) [l2[, as shown in 
Fig. [Id). 

It was shown that the most sensitive states to the dif- 
fusion are the symmetrical CSS with 80 — tt/2, corre- 
sponding to equal populations between the two internal 
states [7| . The phase-diffusion time scale is x^d = 1/ VJ 
(with j = N/2) [2j-|5|. Here, we generalize it for arbitrary 
initial CSS 



|*(0)) = \9 ,<h) = £ C "^'' TO > 
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FIG. 1: Husimi Q function Q(9,<p;i) for initial CSS |0o = 
7r/2, </>o = 0) on Bloch sphere (a), with Bloch vector s = 
(J) = j(sin 6 1 cos <f>, sin # sin 0, cos 9) and j = N/2 = 30. The 
Q functions in the phase space (<j>, s z ) for various times: (b) 
t = 0, (c) t = t s = 0.0764x _1 , and (d) t = 2t d = 0.2651x _1 . 
Time scale x^s — 3 1/ ' 6 (2j) _2 ' /3 is for optimal spin squeezing 
[HQ, and xtd = j~ 1/2 is phase-diffusion time The Q 

function is normalized by its maximum value. 



where the amplitudes 
j - m 



(£)« 



*°(.5) 



In single-particle picture, the CSS corresponds to all the 
atoms occupying a superposed state 15]: cos(#o/2)|2) + 
sin(6* /2)|l), where the polar angle 9 and the az- 
imuth angle <pQ determine atom population and the rela- 
tive phase between the two internal states, respectively. 
In Fig. [2] we plot |c m | 2 as a function of m. One can 
find that the probability distribution of the CSS can be 
treated as a Gaussian wave packet [see below Eq. ([IT])], 
At any time t, the spin system evolves into 



!*(<)) = E 



Croe -i( 5m+xm =)t m) 



(6) 



i.e., a superposition of atomic number state \j, m) = 
\j — m)i\j + m) 2- Due to the presence of atom- atom in- 
teraction (x 7^ 0), each number state has different phase 
evolution rate, which in turn lead to collapse and re- 
vival of the Rabi oscillation, a phenomenon that is well- 
known in quantum optics |14j . Similar effect has been in- 
vestigated in two-component BECs 12411 . atom-molecular 
BECs [H-[27t) and exciton emission 128 301 - 

Phase diffusion of the BECs considered here is in fact 
collapse of the first-order temporal correlation function 

10: 



9® 



\P®\ \(J+)\ 



pM 5 



(7) 



where p$ = (b\bi)/N with k, I = 1, 2 are the elements of 
the single-particle density matrix Q ■ The coherence g^ 



is observable in experiments by extracting the visibility 
of the interference fringes [8 - 12 1 . The expectation value 
(J,) =<*(t)|J,|*(t)), reads 



< = E 



m Cm 



j cos , 



(8) 



which is a constant for a fixed polar angle of the initial 
CSS #o- This is because of conserved operator J z with 
respect to the Hamiltonian H . Similarly, we obtain the 
expectation value (J 2 ) = j 2 cos 2 9q + (j/2) sin 2 9 , and 



— e i(4>o+St+Xt) 



cot E (j : - m ) 



1 2 e 2im X t 



= j e i O°+ <5t ) S i n # Q [cosxt + icos^sinxi] 2 ^ 1 (9) 

where we have used the relation: c rn+ \ = (j — m) 1 / 2 ^ + 
to + l)- 1 / 2 cot(6»o/2)e"^°c m . Inserting Eq. © and Eq. 
(|9]) into Eq. J7|), we further obtain the exact solution of 
the coherence 



5i2 } W= [1-sin 2 (e )sin 2 (#) 



(10) 



which shows a decay of the coherence, i.e., phase diffu- 
sion. Such a dephasing process depends sensitively on 
the self-interaction strength x and the polar angle 9q of 
the initial CSS; while the detuning 5 and the azimuth 
angle O gives vanishing contribution to the coherence. 
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FIG. 2: Probability distribution, |c m J 2 as a function of m for 
the initial CSS \9 ,<f>o}- (a) 9 = tt/2 and (b) 9 = tt/3. The 
empty circle is given by Eq. ((5|, and solid line is given by 
Eq. (pi]). The parameter j = N/2 = 30. 



Time scale of the diffusion is determined by short-time 
behavior of the coherence. Following Imamoglu et al. 0], 
we treat the initial CSS as a Guassian wave-packet with 
its peak located at to ~ ( J z ) = j cos#o> namely 



1 



■ exp 



2(AJ Z ) 2 



[2tt(AJ z ) 2 ]i/2 
where the width 

(AJ z ) 2 = {J 2 z )-{J z ) 2 = (j/2) S m 2 9 



(11) 



(12) 



Due to conserved particle number TV (= 2j), h\ = j — 
J z and hi = j + J z , and t hus atom number variances 
(An x ) 2 



(An 2 ) 2 = (AJ Z ) 2 [17| 
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In Fig. [2] we check the validity of Eq. (fTTI) by com- 
paring it with the exact result of |c m | 2 [where c m is 
given by Eq. (J5J], and find both results fit with each 
other. For symmetric BEC with O = 7r/2, we have 
(J z ) = and atom number variance (A J z ) 2 — j/2; while 
for 0q ^ 7r/2, the variance (width) becomes narrow, as 
shown in Fig. [2jb). Inserting Eq. (fTTj) into Eq. ([9]), and 
replacing the discrete sum over m by an integral, we get 



(J, 



x exp|— p 
where we have set 

9 : 



J_ e i(<Po+St+ x t) e 2i(j z ) x t t ^0 

x / \j - (J z ) - V2(AJ z )p] 

"'Pmin 

.,2 



2iV2 x t(AJ z )p}dp, (13) 



(Jz) 



V2(AJ Z ) 



(14) 



For # ~ tt/2 and j — > oo (large N limit), the integral up- 
per limit p max = Vjtan(#o/2) — > oo and the lower limit 
Pmin = — \/j cot(#o/2) — > — oo. In the short-time limit 
( X t < 1), 1 - 2i X tcos 2 {6 /2) ~ exp{-2^cos 2 (6» /2)}, 
and thus from Eq. (|13[) we obtain 

(J+) ~ j sin(0 o ) exp{i [<j) Q + St + X t(2j - 1) cos O ]} 
xexp [-2 x 2 (AJ 2 ) 2 i 2 ] , (15) 

where two integrals: e~ p e lXp dp = v / 7re _A / 4 and 

pe~ p e lXp dp = iAy / 7re _A / 4 /2 have been used in de- 
rive of the above result. Thus, the short-time solution of 
the coherence reads 



exp 



-2 X 2 (AJ Z )V 



exp 



-{t/Uf , (16) 



with characteristic time scale of the phase diffusion 
1 1 



Xtd 



V2(AJ Z ) v^sin^o 



(17) 



Obviously, the phase-diffusion time X t d = for the 

initial CSS with O = ir/2 @, 0,111; while for 9 ^ tt/2, 
our results show that the diffusion becomes slowly due to 
an enhanced factor sin 0q in the phase-diffusion rate [see 
also Fig. [3]. The revival of the coherence occurring later 
at a time X t r = ir 0, 0] is observable in real experiment, 
such as t r = 7r/x ~ 108.7 microsecond for nonlinearity 
X ~ 2ir x 4.6 Hz and particle number N — 60 [l2| . 
Finally, we note that similar results of Eq. ^ and 
have been obtained in Ref. 31] . However, the 



Eq. {T5J) 

authors focused on the increase of the diffusion time as 
X — > 0, but not 0q ^ tt/2. In addition, we emphasis 
that for negative X the so-called phase separation may 
take place due to dynamically unstable of the BECs 32- 
|39| . Both the phase diffusion and the phase separation 
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FIG. 3: Time evolution of the first-order temporal coher- 
ence (i) for various polar angles: 6o = 7r/6 (dotted line, 
squares), 7r/4 (dashed line, triangles), and 7r/2 (solid line, cir- 
cles). The empty squares, triangles, and circles are given by 
Eq. ()10[) for corresponding 8q's, while the curves are given 
by Eq. (|16p . Horizontal grid line gives the value 1/e, verti- 
cal lines denote the diffusion times td for different 8's. Other 
parameters are N = 60, \ = 2 7r x 4.6 Hz [T^ ]. 



reduce the Ramsey fringes' visibility (i.e., the first-order 
coherence g^). However, the latter effect results from 
dynamics of spatial degree-of-freedom, rather than that 
of the internal states considered here. 

In summary, we have investigated phase diffusion of 
two-component BEC for arbitrary initial CSS. We show 
analytically that for the CSS with Oq ^ tt/2, the diffusion 
process is suppressed due to the reduced atom number 
variance (AJ Z ) 2 below the standard quantum limit j/2 
[see also Eqs. (TT21 and (|17p]. Our analytic results are 
based upon the SMA, which works well to describe the 
condensates tightly confined in a three-dimensional (3D) 
trapping potential. For a lower dimensional cases [40|, 
the coherence [or the fringe visibility] is expected 
to decay more quickly due to the whole continuum of 
excitation modes in confined degrees of freedom 12J, |41 1 . 
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